AD- A 198  795 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  FACE  (Whmn  Dmta  ISntmrod) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


REPORT  NUMBER 


2.  GOVT  ACCESSION  NOJ 


NA 


2  Ado  Z.iQ-CH _ 

TITLE  (and  Subtftfm) 

Theoretical  Studies  of  the  Chemical  Dynamics  of 
Unimolecular  Processes 


3.  RECIPIENT'S  CATALOG  NUMBER 

NA  y  v-  ,  .  . 


5.  TYPE  OF  REPOlW  A  PERIOD  COVERED 

Final  Report 

June  L,  1985-  May  31,  1988 


S.  PERFORMING  ORG.  REPORT  NUMBER 


»UTHORC«J 

Donald  L.  Thompson 


8.  CONTRACT  OR  GRANT  NUMBERS) 

DAAG29-85-K-0039 


performing  organization  name  ano  aooress 

Department  of  Chemistry 
Oklahoma  State  University 
Stillwater,  OK  74078 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  &  WORK  UNIT  NUMBERS 


NA 


controlling  office  name  ano  address 

J.  S.  Army  Research  Office 
3.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


12.  REPORT  DATE 

July  19,  1988 


<3.  NUMBER  OF  PAGES 

18 


14-  MONITORING  AGENCY  NAME  A  AOORESSfl/  dftfmrmnl  from  Controlling  Oil! cm) 


IS.  SECURITY  CLASS,  (of  Mg  rmport) 

UNCLASSIFIED 


15a.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


IS.  DISTRIBUTION  STATEMENT  (of  Mm  Rmport) 

Approved  for  public  release;  distribution  unlimited. 


DUG 


17.  DISTRIBUTION  STATEMENT  (of  Ola  mbmtrmct  mntmrmd  In  Block  20,  II  dlttmrmnf  from  Rmport > 


NA 


AUG  1  5  1988  n 


H 


18.  SUPPLEMENTARY  NOTES 

The  view,  opinions  and/or  findings  contained  in  this  report  are  those  of  the 
author(s)  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  documentation. 


IS.  KEY  WORDS  (Conflmto  on  tmmmrmm  mldm  it  nmcmmmmrf  tad  Idontifr  by  block  number) 


2a  ABSTRACT  (Tam Raw  om  rmrmrmm  X*  ft  naaaaaacy  mad  tdmntttr  by  block  number) 

\\ 

J  The  purpose  of  this  research  program  was  to  develop  computational  methods  and 
carry  out  calculations  to  investigate  the  chemical  dynamics  of  large  molecules  of  interest 
as  energetic  materials.  Classical  dynamics  studies  were  made  of  the  intramolecular 
vibrational  energy  redistribution  and  unimolecular  reactions  in  a  wide  variety  of 
polyatomic  molecules,  including  benzene,  toluene,  dimethylnitramine,  hydrogen  peroxide, 
HONO,  and  methyl  nitrite.  These  studies  are  described  in  this  report.  ^  ^ 

'continued 


DO  , ' 


JAM  71 


1473  EDITION  OF  I  NOV  SB  IS  OBSOLETE  10 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fNRwt  Dmtm  Bntmrmd) 


SCCU WTY  CLASH  FICATIOH  QW  THIS  PAqtfirh—  Data  Snlmrmif) 


20.  ABSTRACT  CONTINUED 

The  technical  personnel  who  contributed  to  the  research  on  this  project  include 
the  PI,  three  postdoctoral  research  associaties,  seven  graduate  students,  one 
undergraduate  student,  and  a  visiting  student  from  Europe.  At  least  twenty  publications 
in  refereed  journals,  one  MS.  thesis,  and  one  Ph.D.  thesis  resulted  from  the  research. 

The  PI  gave  25  invited  talks  at  conferences,  DoD  laboratories,  and  universities  on  the 
research. 


I 

$ 

8 

$ 

1 


i 


i 

a 


§ 

i 


JVM 


FINAL  REPORT 


Theoretical  Studies  of  the  Chemical  Dynamics  of  Unimolecular  Processes 


AUTHOR  OF  REPORT: 

Donald  L.  Thompson 
Department  of  Chemistry 
Oklahoma  State  University 
Stillwater,  Oklahoma  74078 


PERIOD  COVERED  BY  REPORT:  June  1,  1985  -  May  31,  1988 
ARO  PROPOSAL  NUMBER:  22238-CH 
CONTRACT  OR  GRANT  NUMBER:  DAAG29-85-K-0039 


*  -  *  V  V 


w> 


am; 


TABLE  OF  CONTENTS 


Abstract 


Personnel 


Final  Technical  Report  . 3 


References 


Publications,  Theses,  and  Oral  Presentations  . 16 


Abstract 


SS 

: 

1} 

* 


* 

s 

*•2 

III? 


I 

vl 


s 

§ 

f 


'*S 

I 

s 
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Final  Technical  Report 

We  began  this  research  project  with  the  following  long-term  plan  of  approach: 

First,  detailed  calculations  would  be  performed  for  several  selected  polyatomic  molecules 
in  order  to  further  develop  the  computational  methods  and  obtain  the  understanding  of  the 
dynamics  necessary  to  undertake  reliable,  realistic  studies  of  reactions  in  large 
polyatomic  molecules.  Next,  we  would  carry  out  explicit  dynamics  calculations  for 
reactions  in  molecules  such  as  RDX,  TNT,  and  nitroalkanes. 

In  the  this  grant  period,  we  have  investigated  intramolecular  vibrational  energy 
redistribution  (IVR),  the  fundamental  behavior  of  the  dynamics,  and  unimolecular  reactions 
in  a  variety  of  polyatomics.  We  have  also  explored  various  approximations  and  techniques 
for  improving  the  methods  and  facilitating  the  calculations  for  reactions  in  large  molecules. 
Because  of  the  importance  of  the  underlying  IVR  in  unimolecular  reactions,  it  was  studied 
in  detail.  We  have  made  extensive  studies  of  the  rates  and  pathways  of  energy  transfer 
and  of  unimolecular  reactions  in  several  molecules.  We  have  studied  benzene,  hydrogen 
peroxide,  nitrous  acid,  methyl  isocyanide,  methyl  nitrite,  and  dimethylnitramine.  We  have 
investigated  the  influence  of  molecular  structure,  rotation,  potential-energy  surface 
(harmonic/anharmonic,  couplings,  stretch-bend  attenuation,  etc.),  and  initial  energy 
distributions  on  IVR  and  reactions.  The  productivity  of  the  program  is  summarized  in 
terms  of  publications,  theses,  and  presentations  at  the  end  of  this  report. 

We  have  reported  a  thorough  study  of  intramolecular  energy  transfer  in  DMNA.1 
The  intramolecular  energy  flow  out  local  CH^  stretching  modes  of  DMNA  was  studied  by 
using  classical  trajectories  on  various  potential-energy  surfaces.  These  surfaces  were 
constructed  by  using  the  available  experimental  data  augmented  by  some  limited  SCF-MO 
calculations  using  GUASSIAN-82.  One  of  the  purposes  of  this  study  was  to  determine  the 
importance  of  various  features  of  the  potential;  for  example,  we  wished  to  determine  the 
validity  of  certain  approximations,  such  as  harmonic  functions  for  stretches,  which  might 
be  used  to  reduce  computer  time.  And,  of  course,  the  rates  and  mechanisms  of  the  energy 
transfer  in  DMNA  is  of  basic  interest  and  can  help  in  reaching  a  better  understanding  of 
the  unimolecular  reactions.  Various  initial  conditions  and  three  different  potential- 
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energy  surfaces  were  investigated  in  this  study.  The  difference  in  the  two  anharmonic 
potentials  used  in  the  study  is  in  the  way  the  torsional  potential  was  represented.  In 
one,  the  torsional  forces  were  represented  by  a  truncated  cosine  series  and  in  the  other 
by  nonbonding  potentials  (Lennard-Jones).  In  the  third  potential,  all  the  terms  were 
represented  by  harmonic  functions.  The  potentials  all  give  approximately  the 
experimentally  determined  normal  mode  frequencies. 

The  major  features  of  the  energy  transfer  in  DMNA  are  the  same  for  all  three 
potential-energy  surfaces.  The  energy  flow  is  rapid  and  irreversible  and  occurs  on  a 
timescale  of  approximately  0.25  ps.  The  energy  flow  out  of  the  excited  methyl  stretch 
modes  is  initially  into  the  bending  modes  of  the  same  methyl  group  and  then  into  the  other 
methyl  group,  flowing  first  into  the  bending  modes  and  ultimately  into  the  stretches  of 
the  initially  unexcited  methyl  group.  While  the  energy  transfer  is  not  strongly  dependent 
on  the  potential  energy  surface,  there  are  some  minor  differences  in  the  results  for  the 
various  surfaces.  The  energy  transfer  is  somewhat  slower  for  the  anharmonic  potential  in 
which  the  torsional  forces  are  represented  by  Lennard-Jones  potentials  between  the 
nonbonded  atoms  than  in  the  anharmonic  potential  with  the  torsional  forces  given  by  the 
truncated  cosine  series.  This  comparison  was  made  to  check  the  possibility  that  the 
simpler  potential  (i.e.,  Lennard-Jones)  would  suffice  in  general  studies  of  the  dynamics 
of  DMNA.  Use  of  the  Lennard-Jones  potential  reduces  the  computational  time,  but  is 
probably  less  accurate  in  general.  Thus,  we  have  concluded  that  the  cosine  series  is  the 
more  appropriate  form  to  use.  The  anharmonic  potential  referred  to  in  the  following 
discuss  is  the  one  with  the  truncated  cosine  series  for  the  torsional  potential. 

There  are  some  subtle  differences  in  the  energy  transfer  for  the  anharmonic 

potential  and  the  totally  harmonic  potential.  The  rate  of  flow  out  excited  methyl  stretch 

modes  is  essentially  the  same  for  both  potentials.  There  are,  however,  some  significant 

differences  in  the  details  of  the  energy  redistribution  even  though  the  final  distribution 

of  the  energy  after  0.5  ps  is  about  the  same  for  both  potentials.  Thus,  the  totally 

harmonic  approximation  is  acceptable  if  interest  is  focused  only  on  the  overall  rate  of 

the  loss  of  energy  from  an  excited  methyl  group.  This  conclusion  is  interesting  in  light 

2 

of  the  results  of  our  studies  of  benzene,  in  which  we  found  that  it  is  necessary  to 
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treat  at  least  the  excited  stretch  mode  anharmonically  when  the  potential  function  does 
not  include  interaction  terms  (stretch-stretch  and  stretch-bend).  Such  terms  are  not 
included  in  the  DMNA  potential.  An  even  better  approximation  for  benzene  is  to  treat  both 
the  stretch  and  the  associated  bending  motion  anharmonically.  In  DMNA,  the  totally 
harmonic  potential  (without  interaction  terms')  gives  about  the  same  rate  for  the  initial 
flow  of  energy  from  the  stretch,  but  does  not  give  the  same  results  for  the  flow  of  energy 
out  of  methyl  group  into  the  rest  of  the  molecule. 

The  picture  that  has  emerged  from  this  study  is  similar  to  those  for  other 
molecules.  The  bending  degrees-of-freedom  plan  an  important  role  in  the  intramolecular 
energy  redistribution.  The  energy  flow  out  of  CH  stretch  modes  is  directed  along  a 
specific  pathway  into  the  other  methyl  group.  Rapid  energy  transfer  out  of  an  excited  CH 
stretch  appears  to  be  independent  of  the  molecular  environment. 

The  unimolecular  reactions  of  DMNA  were  also  investigated  by  using  classical 
trajectories  on  a  realistic  potential  energy  surface.^  All  degrees-of-freedom  were 
treated  in  these  calculations.  Values  of  the  rate  coefficient  for  the  simple  N-N  bond 
rupture  reaction  to  give  NO^  were  calculated  from  microcanonical  ensembles  of 
trajectories. 

The  intramolecular  energy  transfer  was  examined  in  the  reactive  trajectories.  The 
energy  flow  into  the  nitro  group  increases  and  rapidly  reaches  a  stable  value  just  before 
the  N-N  bond  is  broken.  The  NO^  bend  initially  absorbs  much  of  the  energy  that  flows  into 
the  nitro  group.  After  the  N-N  bond  begins  to  break  the  energy  in  the  nitro  group  rapidly 
exchanges  between  the  bending  and  stretching  modes  of  the  NC^. 

The  predominant  path  for  unimolecular  reaction  in  the  N-N  bond  rupture.  The  CH 
bond  dissociation  is  precluded  by  the  rapid  flow  of  energy  out  of  the  methyl  group.  There 


is  evidence  that  the  HONO  elimination  reaction  channel  may  play  a  small  role  in  the 
overall  decomposition  of  DMNA  on  our  potential-energy  surface.  Two  percent  of  the 
trajectories  produced  geometric  conditions  in  which  a  HONO  molecule  could  be  formed.  This 
potential  does  not  include  the  HONO  elimination  channel.  Thus,  a  more  general 
potential-energy  surface  that  permits  the  HONO  elimination  channel  was  formulated. 
Preliminary  results  of  the  decomposition  of  DMNA  on  a  potential-energy  surface  that 


includes  the  NN,  CH,  ON,  and  NC  bond  dissociations  and  the  HONO  elimination  reaction  have 

been  calculated.  These  results  indicate  that  the  HONO  elimination  reaction  is  not  very 

important  for  a  random  distribution  of  initial  energy  but  the  rate  of  this  reaction  is 

significantly  increased  by  overtone  excitation  of  one  of  the  CH  stretches.  This  is  an 

interesting  results  and  will  be  the  subject  of  future  study.  The  final  energy  in  the  HONO 

elimination  is  primarily  partitioned  into  vibrational  energy. 

The  isomerization  of  methyl  isocyanide  to  methyl  cyanide  is  one  of  the  more 

thoroughly  studied  unimolecular  reactions.  There  are  several  reason  why  have  recently 
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restudied  it  using  classical  trajectories.  First,  the  intramolecular  dynamics  of 

isomerizations  involving  large  amplitude  bending  motions  is  an  important  basic  problem, 

and  this  is  one  of  the  best  known  examples.  This  reaction  is  particularly  interesting 

since  there  are  CH  bonds  adjacent  to  the  isomerization  reaction  coordinate  which  provide 

sites  for  overtone  excitation.  The  possibility  of  mode  selective  chemistry  has  provided  a 

strong  stimulus  for  much  of  the  recent  research  into  intramolecular  energy  transfer. 

Unimolecular  reaction  rates  and  product  energy  distributions  in  keeping  with  statistical 

assumptions  are  most  common.  However,  theoretical  calculations  have  shown  that  energy 

flow  from  a  mode  often  occurs  along  specific  paths.  This  suggests  that  initial 

nonstatistical  flow  of  energy  may  ultimately  be  used  to  direct  the  course  of  reactions. 

Thus,  we  investigated  this  behavior  in  methyl  isocyanide. 

There  have  been  numerous  experimental  and  theoretical  studies  of  the  isomerization 

of  methyl  isocyanide.  Of  most  interest  to  us  was  the  study  of  the  laser-induced 

isomerization  of  methyl  isocyanide  reported  by  Reddy  and  Berry.  They  used 

laser  spectroscopic  techniques  to  measure  the  isomerization  rate.  They  photoisomerized 

CH^NC  by  directly  pumping  CH  overtone  states.  They  obtained  Stern-Volmer  plots  that  are 

linear,  except  for  the  one  at  v=6.  At  v=6  the  rate  constants  are  a  factor  of  five  times 

larger  than  those  predicted  by  RRK.M  calculations  based  on  thermal  data.  The  broad 

overtone  and  combination  bands  observed  in  CH^NC  were  attributed  to  intramolecular 
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dephasing  and  redistribution  occurring  on  a  timescale  of  less  than  0.1  ps. 

The  isomerization  of  methyl  isocyanide  to  methyl  cyanide  has  been  a  prime  subject 

8 

in  the  testing  of  the  basic  validity  of  RRKM  Rabinovitch  and  coworkers  have  carried 


out  extensive  investigations  of  the  thermal  reaction  of  methyl  isocyanide.  The  thermal 
isomerization  of  methyl  isocyanide  is  well  described  by  RRKM  theory  and  is  considered  a 
prime  example  of  the  success  of  the  theory.  We  have  studied  it  because  of  the  similar  to 
some  of  the  systems  of  interest  in  energetic  materials  and  because  of  the  availability  of 
other  results  with  which  to  compare.  It  is  a  good  case  for  checking  potential-energy 
surface  formulations.  Nitromethane,  and  perhaps  nitrobenzene,  undergo  similar 
isomerizations  to  the  nitrite  form. 

9 

Harris  and  Bunker  carried  out  the  first  trajectory  study  of  the  isomerization 

of  methyl  isocyanide.  Harris  and  Bunker  concluded  that  the  reaction  is  non-RRKM  because  of 

slow  intramolecular  energy  transfer.  The  rate  of  energy  randomization  is  not  competitive 

with  the  rate  of  isomerization.  Bunker  and  Hase1^  also  carried  out  a  detailed 

investigation  of  the  isomerization  of  CH^NC  using  classical  trajectories.  They  concluded 

that  methyl  isocyanide  isomerizes  with  lifetime  distributions  that  are  nonrandom  and 

strongly  dependent  upon  the  initial  conditions.  The  work  by  Bunker,  Harris,  and 
9-10 

Hase  suggest  that  intramolecular  vibrational  energy  redistribution  is  incomplete  on 
the  timescale  of  the  isomerization. 
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Because  of  this  background,  and  the  availability  of  some  new  ab  initio  results  on 
the  potential,  we  re-examined  the  dynamics  of  methyl  isocyanide.  The  overtone-induced 
isomerization  of  methyl  isocyanide  to  methyl  cyanide  was  studied  by  using  classical 
trajectories  on  several  potential-energy  surfaces.  We  have  used  this  system  to  explore  the 
sensitivity  of  unimolecular  reactions  to  various  features  of  the  potential.  The  several 
potential-energy  surfaces  are  variations  of  a  potential  which  we  have  developed  based  on 
available  experimental  and  ab  initio  results. 

The  trajectory  results  for  simple  potentials  which  neglect  stretch-bend 
interactions  show  that  overtone  excitation  of  a  CH  stretch  to  the  v=6  level  at  total 
energies  (including  overtone  excitation  energy)  of  75,  125,  and  150  kcal/mol  does  not 
enhance  the  rate  of  isomerization.  However,  at  an  initial  total  energy  of  200  kcal/mol, 
the  isomerization  rate  is  enhanced  by  as  much  as  a  factor  of  three  by  selective  excitation 
of  a  CH  stretching  overtone.  However,  the  mode  specificity  is  sensitive  to  the 
potential-energy  surface.  When  a  more  realistic  potential  is  used  in  which  the  bending 
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forcc  constants  are  attenuated  as  a  function  of  the  bond  lengths  or  in  which  nondiagonal 
quadratic  coupling  terms  are  included,  the  dominant  reaction  is  CH  bond  dissociation.  The 
rate  coefficient  for  the  CH  bond  dissociation  is  an  order  of  magnitude  greater  than  the 
rate  of  isomerization  at  200  kcal/mol.  The  initial  energy  flow  out  of  an  excited  CH 
stretch  is  rapid  (occurring  on  a  timescale  of  less  than  0.5  ps)  and  is  primarily  into  the 
methyl  bending  modes.  The  energy  that  flows  into  the  bending  modes  does  not  transfer  out 
over  the  time  period  of  5.4  ps  that  the  trajectories  were  followed.  The  methyl  bending 
modes  which  act  as  an  energy  "sink". 

Hydrogen  peroxide  is  a  good  example  of  a  large  polyatomic  molecule,  yet  it  is 

small  enough  to  allow  accurate  theoretical  studies  to  be  done.  Thus,  it  has  been  studied 

extensively  by  several  groups.  We1*  carried  out  a  series  of  studies  of  it  by  using 

various  restrictions  on  the  dynamics  to  determine  behavior  of  overtone  relaxation  and 

reaction.  We  have  also  investigated  the  basic  behavior  of  the  dynamics  for  various 

energies,  geometry  restrictions,  and  restrictions  on  the  motion.  Included  in  these 

studies  are  investigations  of  three  other  systems:  HCCC,  HNNH,  and  HCCH.  In  the  first 

study, 11  we  investigated  relaxation  of  OH  overtones  as  a  function  of  excitation  level 

for  collinear,  planar,  and  nonplanar  models  of  hydrogen  peroxide.  The  relaxation  occurs 

12 

very  slowly  compared  with  that  of  a  CH  stretch  overtone.  The  second  study  was  a 
thorough  investigation  of  the  dynamics  of  hydrogen  peroxide  and  the  other  systems.  Power 
spectra  of  auto-  and  cross-correlation  functions  were  used  to  classify  the  type  of  motion 
for  excited  local  modes  and  to  investigate  the  energy  transfer  mechanisms  at  various 
excitation  levels.  Coherency  spectra  were  used  to  determine  the  mode-mode  interactions  and 
thus  the  energy  transfer  mechanism.  The  onset  of  chaotic  motion  in  hydrogen  peroxide 
occurs  at  about  90%  of  the  energy  required  for  dissociation  of  the  OH  bond,  i.e.,  well 
above  the  limit  for  0-0  bond  rupture.  By  studying  restricted  models,  we  determined  the 
influence  of  various  types  of  modes  on  the  dynamics  and  energy  transfer.  For  example,  it 
was  found  in  general  for  the  four  systems  that  the  chaotic  limit  is  lowered  and  the  rare  of 
energy  flow  from  an  excited  local  stretch  mode  is  enhanced  by  the  presence  of  bending 
degrees  of  freedom.  These  calculations  show  that  the  amount  of  energy  transferred  in  a 
given  long  time  (2  or  3  ps)  from  a  local  mode  is  invariant  to  the  number  of  degrees  of 
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freedom  in  the  chaotic  region  but  it  is  strong  dependent  on  the  number  in  the 
quasiperiodic  region.  However,  the  initial  (for  times  less  than  a  ps)  rate  of  energy 
transfer  in  the  chaotic  region  is  substantially  increased  by  the  number  of  degrees  of 
freedom.  That  is,  there  is  a  difference  in  the  behavior  of  the  energy  transfer  for  the  two 
types  of  motion  on  a  very  short  time  scale,  but  the  energy  transfer  on  the  longer  time 
scales  is  about  the  same. 

We  also  carried  out  a  series  of  calculations  on  unimoleculear  dissociation 

14  15 

hydrogen  peroxide  in  collaboration  with  Dr.  Uzer  at  Georgia  Tech.  ’  In  these  studies 
we  explored  mode  specificity  for  overtone  excitations  of  combinations  of  modes.  We  found 
drastically  different  lifetimes  for  different  combinations  of  modes. 

Recently,  we  have  studied  the  influence  of  rotation  on  the  vibrational  energy 

redistribution  in  hydrogen  peroxide  and  on  the  rate  of  0-0  bond  rupture  for  initial 

conditions  in  which  a  zeroth-order  OH  stretch  overtone  state  is  excited.^  We  find  that 

there  is  a  substantial  effect.  The  rates  of  energy  transfer  and  dissociation  are 

increased  by  factors  on  the  order  of  2.5  by  including  rotation.  We  have  also  recently 

completed  a  study  of  the  effects  of  rotation  on  energy  transfer  in  Benzene*^  and  the 
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HONO  cis-trans  isomerization.  ’  The  effects  of  rotation  on  intramolecular  dynamics 
and  unimolecular  reaction  rates  have  long  been  a  topic  of  considerable  interest.  For 
example,  a  number  of  papers  have  been  written  on  it  with  reference  to  RRKM. 

HONO  is  an  important  product  in  the  initial  decomposition  of  some  molecules  of 
interest  as  energetic  materials.  It  is  also  interesting  from  a  basic  research  point  of 
view.  The  results  for  hydrogen  peroxide  have  provided  a  lot  of  insight  into  the 
intramolecular  dynamics,  energy  redistribution,  and  dissociation  Thus,  it  seemed 
reasonable  to  extend  these  studies  to  HONO.  It  is  one  of  the  simplest  molecules  that  can 
undergo  cis-trans  isomerization  and  is  more  general  than  HOOH  since  it  is  not  symmetric. 
The  barrier  to  isomerization  is  about  9.7  kcal/mol.  Thus,  it  provides  a  real  system  for 
which  questions  of  mode  specific  excitation  influence  on  reaction  can  be  explored.  Based 
on  the  results  for  HOOH  and  experimental  studies  on  similar  systems  by  Bauer  and 


coworkers  ,  it  seems  reasonable  that  the  intramolecular  conversion  in  HONO  would 
provide  a  realistic  situtation  in  which  nonstatistical  effects  might  be  important. 


Pimentel  and  coworkers,^  and  McDonald  and  Shirk^  have  also  published  some 

experimental  studies  that  suggest  that  HONO  will  display  interesting  dynamics  effects. 

22 

Darsey  and  Thompson  carried  out  ab  initio  SCF-MO  calculations  to  obtain  an  accurate 

torsional  potential  function  for  HONO.  Completely  optimized  geometries  using  the  4-31G 

basis  set  were  used  to  calculate  a  potential  which  was  least-squares  fit  with  a  cosine 

series  for  use  in  dynamics  calculations.  The  barrier  to  isomerization  on  the  computed 

potential  is  9.68  kcal/mol.  Classical  trajectories  were  then  used  to  investigate  cis-trans 

17 

isomerization  and  intramolecular  vibrational  energy  redistribution.  The  influence  of 

various  initial  normal  mode  excitations  were  studied.  Power  spectra  were  used  to  examine 

the  dynamics.  The  rate  of  cis-trans  conversion  is  strongly  dependent  on  the  site  of 

initial  excitation.  The  rates  for  cis-trans  isomerization  is  significantly  higher  than 

those  for  trans-cis.  The  results  show  that  the  OH  stretching  mode  is  only  weakly  coupled 

to  the  other  modes.  Energy  transfer  out  of  overtones  of  the  OH  stretch  is  slow  and 

preferential  partitioning  of  initial  energy  in  that  mode  leads  to  a  reduced  rate  of 

isomerization  compared  with  other  initial  energy  distributions  at  the  same  total  energy. 

The  HON  bend  mode  interacts  with  the  torsion  more  strongly  than  do  any  of  the  other  modes. 

Recently,  we  have  completed  a  study  of  the  influence  of  rotation  on  HONO 
1 8 

isomerization.  By  using  a  potential  with  no  coupling  terms,  we  can  show  that  there 
are  strong  effects  of  the  rotation  on  the  coupling  of  the  kinetic  energy.  This  work  is 
being  prepared  for  publication. 

The  intramolecular  vibrational  relaxation  and  the  high-energy  spectroscopy  have 
been  studied  more  thoroughly  for  benzene  than  most  other  molecules.  From  a  theoretical 
standpoint,  it  has  served  as  prototypical  system  for  investigating  the  rates  and 
mechanisms  by  which  energy  initially  localized  in  CH  stretch  mode  becomes  redistributed. 

A  number  of  people  have  reported  theoretical  studies  of  the  IVR  in  benzene.  Because  it  is 
well-studied  by  various  methods,  it  is  an  excellent  molecule  to  use  to  investigate  various 
aspect  of  computational  treatments  of  the  full  dynamics  of  large  molecules.  We  have 
carried  out  a  number  of  studies  of  benzene  aimed  primarily  at  exploring  the  details  and 
methods  of  classical  dynamics  treatments  of  large  molecules.  For  example,  we  have  used  it 
to  check  various  approximations  for  potentials. 


■v.’vv*  uku*  vrnif1  ir^trTjfT^ 


A  variety  of  mechanisms  have  been  advanced  to  explain  the  rapid  rate  of 

intramolecular  vibrational  relaxation  observed  subsequent  to  overtone  excitation  of  CH 

bonds  in  benzene.  Resonance  interactions,  especially  the  Fermi  resonance  between  an 

excited  CH  bond  and  the  associated  CCH  wag,  have  generally  been  advanced  as  instrumental 

23 

in  enchancing  the  rate  of  energy  transfer  from  the  excited  CH  bond.  Shi  and  Miller  and 

24 

Cerjan,  Shi,  and  Miller  have  reported  theoretical  studies  of  lineshapes  of  CH  stretch 

overtones  in  benzene  in  which  they  discuss  the  decay  of  energy  from  the  excited  CH  bond  in 

terms  of  resonance  interactions.  Several  other  theoretical  studies  have  also  reported 

25  26 

numerical  techniques  for  calculation  of  CH  overtone  lineshapes  in  benzene.  ’  Sibert 
27 

et  al.  have  reported  classical  and  quantum  mechanical  studies  that  address  the  decay 

of  energy  from  overtone  excitations  of  CH  bonds  in  the  context  of  resonance  interactions. 

28 

Stannard  and  Gelbart  reported  a  theoretical  study  of  intramolecular  vibrational  energy 

29 

redistribution  in  benzene.  Reddy,  Heller,  and  Berry  proposed  and  tested  several 

possible  mechanisms  for  energy  flow  from  overtone  exitations  of  CH  bonds.  Holme  and 
30 

Hutchinson  have  noted  that  in  polyatomic  molecules  intramolecular  vibrational 

relaxation  occurs  much  faster  than  does  excitation.  As  energy  is  absorbed  by  a  local  mode, 

it  is  almost  simultaniously  relaxing  into  eigenstates  which  project  onto  that  local  mode. 

This  is  evident  in  the  classical  trajectory  calculations  by  the  immediate  rapid  flow  of 

31 

energy  from  excited  zeroth-order  states.  Nagy  and  Hase  have  reported  a  classical 

32 

trajectory  study  of  benzene.  More  recently,  Clarke  and  Collins  have  reported 

trajectory  studies  of  benzene  and  perdeuterobenzene. 

We  have  reported  several  classical  trajectory  studies  of  CH  overtone  relaxation  ir 
2  16  33-35 

benzene.  ’  ’  These  studies  focus  on  approximations  in  potential-energy  surface 

formulations.  The  purpose  was  to  learn  more  about  where  approximation  can  be  made  to 

simplify  the  calulation  of  dynamics  of  large  polyatomic  molecules.  In  the  first  two 

2  33  33 

studies,  ’  potentials  without  coupling  terms  were  used.  In  the  first  study,  we 

computed  trajectories  for  planar  benzene  on  two  types  of  potential-energy  surfaces, 

harmonic  and  anharmonic,  by  using  the  analytic  form  of  Nagy  and  Hase,  but  with  different 

2 

values  of  the  parameters.  In  the  second  study,  the  effect  of  anharmonicity  on  the  IVR 
was  more  fully  investigated.  We  investigated  potential-energy  surfaces  with  varying 


1 
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degrees  of  anharmonicity  to  determine  which  m<  des  can  be  treated  harmonically  without 
affecting  the  rate  of  energy  transfer.  The  results  show  that  it  is  necessary  to  treat  at 
least  the  excited  stretch  mode  anharmonically,  and  a  better  results  is  obtained  if  the  wag 
mode  (which  interacts  most  strongly  with  the  stretch)  is  also  treated  anharmonically. 

These  results  have  ramifications  for  reaction  path  hamiltonian  approaches  where  a  reaction 
mode  is  coupled  to  a  harmonic  bath. 

Following  these  studies,  we  investigated  to  accuracy  of  a  harmonic  force  field 
with  coupling  terms  as  well  as  the  influence  of  out-of-plane  motion  on  IVR.  The  third 
study  *  employed  harmonic  potential  of  the  form  of  the  ab  initio  surface  developed  by 
Pulay,  Forgarasi,  and  Boggs.^  The  results  obtained  in  this  study  indicate  that  energy 
redistribution  in  benzene  is  very  rapid,  occurring  on  a  timescale  of  approximately  .2  to 
.3  ps.  There  appears  to  be  no  localization  of  energy  within  the  molecule  on  more  than  a 
subpicosecond  timescale.  This  is  an  interesting  result  because  the  potentials  used  in 
these  calculations  are  harmonic.  In  our  previous  studies  of  benzene  using  simplier 
(primarily  fewer  coupling  terms)  analytical  forms  for  the  potentials  it  was  found  that 
anharmonic  potentials  were  necessary  in  order  to  correctly  predict  the  rapid  energy 
transfer  out  of  excited  modes.  The  results  for  the  Pulay  et  al.  potential  show  that  the 
importance  of  anharmonicity  depends  on  other  features  of  the  potential-energy  surface.  A 
purely  harmonic  potential  can  predict  the  correct  rate  of  energy  transfer  in  benzene. 
Obviously,  the  harmonic  approximation  is  not  accurate  in  some  respects  at  these  levels  of 
excitation.  Thus,  some  of  the  computed  properties  of  the  energy  transfer,  such  as  the 
details  of  the  pathways,  may  not  be  completely  realistic.  Nevertheless,  it  is  important 
to  note  that  the  harmonic  force  field  gives  the  correct  behavior  for  the  initial 
relaxation  of  the  CH  stretch  overtone. 

We  also  investigated  the  energy  transfer  pathways  by  following  the  energy  in  the 
various  normal  modes  as  a  function  of  time.  There  are  two  distinct  sets  of  ring  modes,  the 
essentially  nonparticipating  low  frequency  modes  and  the  higher  frequency  modes  which  are 
integrally  involved  in  the  rapid  decay  of  energy  from  the  CH  bond.  Within  the  higher 
frequency  range,  participating  set  of  modes  it  appears  that  the  modes  which  possess  the 
most  CCH  bend  character  are  more  strongly  coupled  to  the  excited  CH  bond.  Energy  flows 


preferentially  into  these  modes  and  then  decays  into  other  ring  modes  which  possess 

significant  ring  stretch  and  deformation  character. 

The  effect  of  out-of-plane  motion  on  IVR  was  studied  by  comparing  results  for 

planar  and  nonplanar  benzene.  It  was  found  that  the  inclusion  of  out-of-plane  motion  does 

not  appear  to  enhance  the  rate  of  energy  flow  from  an  excited  CH  bond. 

We  also  used  classical  trajectories  to  investigate  energy  relaxation  of 

various  excited  C-C  stretching,  CCH  bending,  and  CCC  bending  normal  modes  in  benzene. 

These  normal  modes  were  initially  excited  up  to  about  10,000  cm  *.  The  calculated 

relaxation  of  the  excited  C-C  ring  breathing  mode  is  in  agreement  with  the  recent 

37 

experimental  results  of  Chernoff,  Myers,  and  Pruett.  In  general,  the  relaxation  of 
the  excited  normal  modes  involving  C-C  stretching,  CCH  bending,  and  CCC  bending  are  at 
least  several  times  slower  than  the  relaxation  of  the  CH  stretch  local  mode.  Mode-to-mode 
energy  transfer  was  also  investigated.  Only  those  modes  whose  frequencies  are  close  to 
that  of  the  excited  mode  participate  in  the  energy  transfer. 

38 

Cis-irans  isomerization  and  IVR  in  methyl  nitrite  has  also  been  studied.  We 

have  computed  an  ab  initio  torsional  potential  and  fit  it  with  a  cosine  series  as  we  did 

39 

in  the  case  of  HONO.  Initial  conditions  for  microcanonical  distributions  of  initial 
energy  and  CH  stretch  overtones  were  investigated.  We  found  that  the  CH  overtone 
excitations  give  the  same  rates  as  does  the  equipartitioning  of  the  initial  energy. 

We  arc  currently  studying  the  competing  simple  bond-rupture  reactions  to  give  CH^O  +  NO 
and  CH3  +  N02- 
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